Measurement of the current-phase relation of superconducting atomic contacts 
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We have probed the current-phase relation of an atomic contact placed with a tunnel junction in a 
small superconducting loop. The measurements are in quantitative agreement with the predictions 
of a resistively shunted SQUID model in which the Josephson coupling of the contact is calculated 
using the independently determined transmissions of its conduction channels. 

PACS numbers: 74.50.+r, 74.45.+C, 74.78.Na, 73.23.-b, 73.63.-b 



The Josephson effect is a striking signature of extended 
quantum coherent states of matter, as found in superflu- 
ids, superconductors, and atomic Bose-Einstein conden- 
sates. It appears when a weak link allows particles to 
flow between two reservoirs of such quantum systems, 
thereby establishing phase coherence between the two 
corresponding macroscopic wave functions. The effect 
was first predicted and observed for the case of super- 
conductors, where DC electrical supercurrents flow in ab- 
sence of any voltage and AC supercurrents appear under 
a constant voltage bias ll| . Since then, it has also been 
explored in superfluids |2J and in Bose-Einstein conden- 
sates [3j. Although in the field of superconductivity a 
large variety of weak links has been used (tunnel junc- 
tions, proximity effect bridgeSjPoint contacts, graphene, 
carbon nanotubes, . . . ) [a, 0, @|, the basic effect is 

generic and a unifying picture, able to treat on the same 
footing all the different coupling structures, emerged in 
the framework of mesoscopic superconductivity. Within 
this framework, the basic Josephson weak link is a sin- 
gle conduction channel of arbitrary transmission proba- 
bility r 0] connecting two superconducting electrodes. 
For a channel shorter than the superconducting coher- 
ence length, the Josephson coupling between both sides 
is established by a single pair of "Andreev bound states" 



[id These are described as resonant electron- hole 
quasiparticles states spreading in both electrodes, with 

energies E±(6,t) = ±A (l - Tsin 2 (<5/2)) 1/2 , A being the 
superconducting gap, and 5 the phase difference between 
the order parameters on both sides. These states carry 
opposite supercurrents 1^(5) — 2-K(j)Q 1 dE^r/dS, where 
4>a = h/2e is the flux quantum. At zero temperature 
only the lower state is occupied, and the current-phase 
relation for the weak link is simply I~ (S) . The criti- 
cal current of a channel 7° = max{/~(<5)} is the max- 
imum supercurrent that can be sustained in absence of 
fluctuations of the phase. Beyond this value, a voltage 
develops across the system, i.e. transport becomes dissi- 
pative, and is perfectly understood in terms of multiple 
Andreev reflection (MAR) processes 12 1, which also de- 
pend strongly on r. In general, any phase coherent con- 
ductor can be described as a collection of independent 
conduction channels, characterized by its set of trans- 



mission coefficients {t^}, and the global current- phase 
relation J/ Ti \(<J) = ^Zjr(^) for an arbitrary weak link 
is simply the sum of the contributions of its channels. 
It is then clearly of fundamental interest to measure the 
current-phase relation for a single channel of arbitrary 
transmission. Atomic contacts of superconducting met- 
als are suitable systems to test these ideas, even if they 
tend to comprise not just one but a few channels [T3|. A 
measurement of the current-phase relation of atomic size 
point contacts was first performed by Koops et ai. [l4j . 



but a quantitative comparison with theory could not be 
performed because the {r^} were not known. Since then, 
a reliable method to determine the transmissions was de- 
veloped [H|], based on the measurement of the dissipa- 
tive MAR current under a voltage bias. In this Letter 
we present measurements of the current-phase relation 
of well characterized aluminum atomic contacts, and a 
direct comparison with theory, with no adjustable pa- 
rameter. 

The principle of our experimental setup is shown 
schematically in Fig. la, and is designed to allow for both 
voltage and phase bias of the samples, giving this way 
independent access to the transmission probabilities and 
to the current-phase relation respectively. The sample 
consists of an atomic contact (phase difference ft) and a 
tunnel junction (phase difference 7) embedded in parallel 
in a small superconducting loop [16[, hence forming an 
asymmetric SQUID, as shown in the scanning electron 
microscope image of Fig. lb. Note that similar ljj or 
related [181] setups are being used in other laboratories. 
The atomic contact is obtained using a microfabricated 
break junction [l9l | , and is characterized by a critical cur- 
rent I? T .\i typically a few tens of nA, much lower than the 
critical current of the tunnel junction Iq w 740 nA. The 
SQUID is placed in parallel with an on-chip rC circuit, 
which dominates the parallel impedance of the line and 
controls the phase dynamics. The sample, which is fab- 
ricated on a metallic substrate coated with a polyimide 
layer, is thermally anchored to the mixing chamber of a 
dilution refrigerator. The bias current It is provided by a 
room temperature voltage source connected through a se- 
ries of 50 fl attenuators placed at low temperatures, and 
a discrete macroscopic resistor R mounted at the same 
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FIG. 1: (color online) (a) Schematic experimental setup: DC 
SQUID formed by an atomic contact (phase S) and a tun- 
nel junction (phase 7). The on-chip capacitor C « 21 pF 
is formed between the metallic substrate and a 100 nm-thick 
gold electrode (1.3mm 2 ), with a 1.6 fjxa- thick dielectric poly- 
imide layer. The resistor r ~ 0.6 f2, corresponds to the sheet 
resistance of the capacitor gold electrode. The bias current 
is governed by voltage source U and discrete macroscopic re- 
sistor R m 25 mounted close to sample, at base temperature 
(20 mK) of dilution refrigerator, (b) SEM image of SQUID 
loop. Tunnel junction is fabricated using double-angle evap- 
oration of aluminum through a suspended mask. This results 
in a parasitic structure of no practical importance. Bright 
regions on left corners correspond to the gold thin films that 
connect the superconducting loop to the rest of the circuit, 
and provide the top plate of the capacitor. These normal 
regions also act as quasiparticle traps. 



temperature as the sample. Voltage and current in the 
SQUID are both measured with low-noise voltage ampli- 
fiers, the latter from the voltage drop across the resistor 
R. 

The idea behind this setup is twofold [la |. On the 
one hand, it allows to obtain the dissipative part of the 
current- voltage characteristic of the contact Ir T .\(V) — 
^SQuro(^) — Ijiy), as the difference between the one of 
the SQUID /squid {V) and the one of the tunnel junction 
Ij{V) alone, as shown in the lower panel of Fig. 2. The 
latter is measured after fully opening the break junction 
[2^ |. One then determines the transmission probabili- 
ties {n} and the gap A ~ 180 ueV, by fitting I{ Ti }(V) 



with MAR theory 15, 2(J 21, 22], as shown in the upper 
panel of Fig. 2 for three different contacts. On the other 
hand, on the supercurrent branch of the SQUID (see in- 
set in Fig. 2), it is possible to impose a phase difference 
on the contact using both the external flux and the cur- 
rent bias as control knobs, and to use the tunnel junction 
as a threshold detector to measure the current flowing in 
the loop 12711 . Indeed, the loop is designed to be small 
enough [28j so that, to a very good approximation, the 
two phases are linked by the magnetic flux <f> threading 



the loop, according to 8 — 7 = ip = 2tt(/)/<Po and therefore 
lb = Iq sin 7 + I{ Ti }{l + </?)■ In the limit Iq S> J? Tf i the 
critical current I c of the SQUID should be reached when 
7 ~ 7r/2, and therefore its variations with the external 
flux are I c (</?) ~ Jo + Is Ti \(ip + tt/2). Therefore, the pe- 
riodic modulations of I c (ip) around the critical current 
of the junction probe directly the current-phase relation 
of the atomic contact. It is however important to note 
that in practice, due to fluctuations, both quantum and 
thermal, the system "switches" stochasticaly from the 
supercurrent branch to the dissipative branch before the 
bias current reaches I c . This switching process is char- 
acterized by a rate T. It is nevertheless still possible, as 
we show hereafter, to probe the current-phase relation of 
the contact from measurements of the switching current 
of the whole device as a function of the magnetic flux. 
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FIG. 2: (color online) Upper panel: Experimental current- 
voltage curves for three atomic contacts (symbols) for < 
eV < 2A, compared to best fits J {t . } (V) using MAR (full 
fines). Fits provide gap A ~ 180 ^teV and set {n} of transmis- 
sion coefficients for each contact: AC1 — > {0.62; 0.22; 0.07}; 
AC2 -» {0.957; 0.19}; AC3 -> {0.993; 0.14}. Lower panel: 
I(V) curve for contact AC3, obtained as the difference be- 
tween Isquid(V) and Ij(V). Inset: Large scale character- 
istic of SQUID, displaying supercurrent branch at V = 0. 

Typically, we apply 10 4 bias current pulses of ampli- 
tude s = Ib/Io and duration r p ~ 40 /is, and measure 
the switching probability P(s) = 1 — e~ r ( s ) Tp as the ra- 
tio between the number of switching events and the total 
number of pulses. For each value of (p we adjust the 
current pulse amplitude s*(p) so as to keep a constant 
switching probability P(s*) = 0.6 (corresponding to a 
rate T* — 23.3 kHz), which leads to the best sensitivity 
with respect to flux response. The s*(p) curves measured 
in this way for the three contacts of Fig. 2 are shown 
as symbols in Fig. 3. As the absolute value of the flux 
through the loop is not known, and the current measure- 
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FIG. 3: (color online) Symbols: measured switching cur- 
rent (s*(<p) — so)Io as a function of applied flux 4>/<j>o, for 
the three SQUIDS corresponding to the contacts of Fig. 2. 
Curves AC3 and AC1 shifted for clarity. Dashed curves: 
predicted ground state current-phase relation 7T_ j(<5). Full 
lines: predictions of resistively shunted SQUID theory at 
Tcsc = 130 mK on the basis of Eq. (1). The transmission sets 
indicated in Fig. 2 caption have been used for both theories. 



junction, with Ej = (f) I Q /2ir, the second one is the en- 
ergy arising from the coupling to the current source, and 
the last term is the Josephson coupling introduced by the 
atomic contact. Figure 4 shows U-(j) for a SQUID with 
a single channel contact (t = 0.99), and for comparison 
the potential of the tunnel junction alone. 




FIG. 4: (color online) Full line: washboard potential of a 
SQUID with a single channel contact for r = 0.99, s = 0.87 
and (p = 0, as function of the Josephson junction phase 7. 
Thermal activation allows the phase to escape at a rate V 
above the barrier of height AU. Dashed line: washboard 
potential of Josephson tunnel junction alone for the same pa- 
rameters. 



ments suffer from small drifting offsets, the experimental 
curves had to be shifted in both directions by an arbitrary 
amount in order to compare them with the theoretical 
predictions. Note however that the vertical shift corre- 
sponds essentially to the average value (s*(ip)} , which 
is very close to the switching current sq of the junction 
alone under the same conditions. Figure 3 also shows the 
calculated current-phase relations I{ Ti }{<p + ir/2) for the 
corresponding sets {r^}. There is an overall qualitative 
agreement between the experimental data and these sim- 
ple predictions. The discrepancies are significant only 
for contacts AC2 and AC3, which both contain a highly 
transmitted channel, and arise mainly around a phase 
5 = 7r. These differences can be understood almost com- 
pletely by taking into account the phase fluctuations im- 
posed by the dissipativc elements of the electromagnetic 
environment in which the SQUID is embedded. It is well 
known that in such a dissipative biasing circuit the phase 
across the SQUID is a dynamical variable governed by a 
Langevin equation, equivalent to the one obeyed by the 
position of a massive particle evolving in a "tilted wash- 
board potential" in presence of friction [2{|. Assuming 
that only the ground Andreev state of each channel of 
the atomic contact is occupied, the total potential of the 
SQUID is given by: 

U-{l) = -Ej cos 7 - Ejs^ - 'YjE- (7 + if,n) (1) 

i 

where the first term is the Josephson energy of the tunnel 



The overall shape of the potentials is qualitatively the 
same but for very highly transmitted channels (r > 
0.999), and the physics is therefore similar to the well 
known case of tunnel junctions. For the actual parame- 
ters of the setup, one can neglect quantum fluctuations 
and treat 7 as a classical variable. For < s < 1, 
the equivalent particle oscillates around a local mini- 
mum of the potential at the plasma frequency lu p (s) — 
lo {1 - s 2 ) 1 / 4 , with u ~ (27r/ /^ C) 1 / 2 [3Cj. The tilt 
of the potential increases with s, and the thermal energy 
ksT becomes eventually comparable to the potential bar- 
rier height At/ (7) = U-(j max ) - t/_(7 mln ), where 7 mi „ 
(7max) is the phase at which the potential presents a lo- 
cal minimum (maximum) . The particle can then be ther- 
mally activated over the barrier and escape from the well 
at a rate 

r( a ,<p) ~ ^ e -AC/ (sw )/fe fl T (2) 

before the current actually reaches the critical current. 
The biasing circuit is such that, once escaped, the parti- 
cle runs away indefinitely and a voltage suddenly devel- 
ops at the edge of the SQUID, according to V = 0o7/27r. 
This corresponds to the "switching" detected in the ex- 
periments. 

We first performed switching measurements of the 
Josephson junction alone, which is a well known case 
[29j ]. The reduced bias current corresponding to the im- 
posed escape rate is of the order of s = 0.87, correspond- 
ing to a phase 7 m i n = 70 = arcsin(so) — 0.67 (tt/2). 
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The s-dependence of the switching rate agrees precisely 
with Eq. (2), and yields the escape temperature T csc w 
125 mK. Although it is significantly higher than the re- 
frigerator temperature T = 20 mK, showing that the 
electrons in the dissipative elements of the biasing cir- 
cuit are heated by some remaining spurious noise, it does 
not hinder the measurement of the current-phase rela- 
tion. To compare the experimental data for the SQUIDs 
with theory, we assume this measured T osc to be the ac- 
tual temperature determining the phase fluctuations in 
all cases, and given the exponential dependence of the 
rate on AJ7(s, ip), that a constant escape rate corresponds 
to a constant barrier height. We then use Eq. (1), with 
the {Ti} obtained independently (Figj2j) , to calculate, for 
each SQUID, the current s*(ip) leading to the imposed 
rate T* . The curves calculated in this manner are shown 
as full lines in Fig. 3 and describe the experimental data 
significantly better than the T — theory. Note that this 
procedure assumes that in each channel only the ground 
Andreev state is occupied, and that the only effect of the 
finite temperature is on the dynamics of the phase. This 
is not forcedly the case for the contacts measured here, 
as fcsT esc is not much smaller than the minimum energy 

1 /2 

gap 2A (1 — r) ' between Andreev levels at 5 = n, but 
the resulting population of the excited state is still too 
small to have a detectable consequence. 

In conclusion, we have measured the current-phase 
relation of superconducting atomic contacts covering a 
wide range of transmission coefficients, and accounted 
quantitatively for the results using the mesoscopic the- 
ory of the Josephson effect. Note that the experiments 
described here probe just the ground Andreev state of 
each channel but that it should be possible and interest- 
ing to also probe the excited state, through microwave 
spectroscopy for instance. Furthermore, by designing a 
proper environment such that the life-time of this excited 
state and the dephasing time between the two states are 
long enough, one could envision to create coherent quan- 
tum superpositions of them [31 1 . 
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